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SmTARY 

This  report  explores  tlaree  approaches  to  the  turning 
point  problem  for  linear  equations. 

The  fir.jt  approach  exhibits  a  full  set  of  solutions 
of  the  equation  as  contour  integrals.   This  is  illustrated 
hy   the  triangular  system  and  the  hyper^eometric-type  equation. 

The  second  approach  regards  the  sr/stem  as  being  a 
perturbation  of  a  comparison  equation.   In  particular,  if 
the  leadi.Mg  matrix  of  the  comparison  equation  is  dia^onalizable , 
the  rapid  convergence  of  the  Picard  iterant s  is  shoirni  by  an 
elaboration  of  a  scheme  due  to  ll.  0.  Friedrichs. 

The  third  approach  regards  the  system  as  bein^  a 
composite  system  resulting  from  the  amalgamation  of  component 
solvable  systems.   Particular  attention  is  given  to  the 
reduction  attainable  due  to  a  non-coalescing  eigenvalue. 
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I .   Introduction 

The  turning  point  problem  has  received  onl"  sparse 
treatment  for  the  general  linear  equation.  However,  some  of 
the  approaches  employed  for  special  equations  appear  to  be 
readily  extensible  to  more  general  situations.   It  is  the 
intent  of  this  report  to  indicate  some  of  the  possibilities 
in  this  direction  which  were  explored  during  the  past  year. 
The  '.;or!-:  is  still  in  its  preliminary''  stage  ::o   that  the 
discussion  is  suggestive  rather  than  definitive. 

The  turning  point  problem  and  the  related  irregular 
singular  point  problem  is  formulated  in  section  II,  This  is 
folloT/ed  by  a  discussion  in  section  III  of  special  equations 
for  TJhich  the  turning  point  problem  is  directly  solvable  by 
contour  integrals.   The  first  such  equation  is  the  often 
neglected  triangular  system  (vriiere  the  equation  for  each 
component  of  a  vector  solution  is  decoupled  from  all  the 
preceding  components).  The  other  solvable  equation  considered 
at  length  is,  in  essence,  the  equation  completely  solvable 
by  the  generalized  hypergeometric  functions.   The  discussion 
of  this  equation  is  based  upon  the  work  of  H.  Scheffe  [1]""' 
and  H.  L.  Turrittin  [1].   Incidentally,  both  of  the  above 
equation  t-jz^es   are  solvable  even  if  the  turning  point  occurs 
at  a  regular  singularity. 


,  This  notation  refers  to  the  first  paper  listed  under  the 
name  of  Scheffe  in  the  bibliography  at  the  end  of  this  report. 


'Jhen  the  equation  is  not  directly  solvable  one  attempts 
a  perturbation  scheme  which  goes  back  to  Liouvllle  and 
Steckloff .   This  scheme  regards  the  differential  equation  as 
being  a  perturbation  of  a  comparison  equation.  The  differential 
equation  is  then  converted  into  a  Volterra  integral  equation 
v;h03e  inhomogeneous  term  is  the  solution  of  the  comparison 
equation.  'Tiat  remains  is  to  demonstrate  the  rapid 
convergence  of  the  iterant  solution  of  the  integral  equation. 
One  such  demonstration  is  due  to  K.    0.  Friedrichs  [1,  p.  191] 
and  this,  denoted  as  the  strong  perturbation  approach,  is 
elaborated  in  section  IV.  Another  demonstration  stems  from 
R.  E.  Langer  [1]  and  a  possible  extension  of  this  weak 
perturbation  is  sketched  in  section  V.  Finally,  one  may 
regard  a  system  as  being  a  composite  system,  resulting  from 
the  amalgamation  of  component  solvable  systems.  This  is 
indicated  in  section  VI  v;ith  particular  attention  to  a  system 
in  which  an  eigenvalue  does  not  coalesce  at  the  turning  point 
(Langer  [2]). 

VJhatever  merit  this  report  may  possess  is  due  to  the 
valuable  information,  advice,  criticism,  and  encouragement 
which  was  unstintingly  given  by  Professor  K.  0.  Friedrichs. 


Notation 

Tlirouchout  the  report  the  follov;lng  notation  1: 
maintained: 


denotes  differentiation  v;lth  respect  to  x 
X   is  the  complex  independent  variable 
X  =  0  is  the  turning  point 

p  =  e"   is  a  large  paraneter  (taken  real  for  convenience) 
A,  or  any  other  capital  letter  in  an  equation  is  a 
square  matrix 


ie=^^V 


II.   The  Turning  Point  Problem 

Consider  the  system  of  first  order  equations  for  the 
matrix  Y 

Y'  =  A(x,e)Y 

i-zhere  e  is  a  parameter.   If  A  is  a  'suitable'  function  of  x 
and  e  then  it  is  well  knov/n  that  as  e  — *  0,  there  e::ists  a 
full  compleraent  of  formal  solutions  Y(x,e)   (Noaillon  [1], 
Tiirrlttin  [2]).   Furthermore,  in  a  certain  'minimum  sector'  of 
the  complex  x-plane,  boimded  by  two  'transition  lines'  which 
meet  at  a  ' turning  point',  there  exists  an  actual  matrix 
solution  v/hich  is  'asymptotic'  to  the  formal  solution 
(Birlchoff  [1],  Birkhoff-Langer  [1],  Trjitzinsky  [2], 
Htikuhara  [1],  V.fasow  [1,  chap.  V]).  Moreover,  this  actual 
solution  can  be  shown  to  preserve  its  asymptotic  form  beyond 
this  'minimum  sector'  up  to  a  certain  'maximum  sector'  which 
includes  several  transition  lines  in  its  interior  (Priedrichs 
[1,  p.  155].   Ifj  as  is  usual,  this  maximum  sector  is  the 
sector  of  largest  extent  for  the  ^iven  asymptotic  form  then 
it  is  denoted  as  the  Stokes  sector  and  its  bounding  lines  are 
called  Stokes  lines  (Stokes  [1]).   Of  course,  in  a  decoupled 
system  the  Stokes  sectors  are  in  general  larger  than  the 
maximiim  sectors.   In  any  event,  the  neighborhood  of  the  turning 
point  can  be  covered  by  a  set  of  overlapping  maximum  sectors, 
say  Jo  .,  in  each  of  v/hich  there  is  an  actual  solution,  say  Y.  > 
such  that 


'Hvj     'it' 


Y^  ^  Y   in 


A  ■ 


The  problem  Is  to  find  the  asymptotic  behavior  of  each 
Y.  in  all  the  other  sectors. 

To  solve  this  problem  it  is  sufficient  to  find  the 
asymptotic  form  (with  a  remainder  the  order  of  a  power  of  e) 
of  each  Y.  at  the  turning  point  v/ith  the  proviso  that  the 
colxiTiins  of  this  asymptotic  form  are  independent. 

One  then  completes  the  argument  by  noting  that  due  to 
the  linearity  of  the  equation,  any  two  fundamental  solutions, 
say  Yu  and  Y,^,  can  be  related  through  the  mediation  of  an 
arbitrary  ordinary  point  x  by  the  equation 

^■^h(^)'^h'(^o)  =  V^^^'-^^o)   ' 
i.e., 

Y^(x)  =  Y^(x)K(e)   where  K( e )  =  YJ^(Xq)Y^(x^) 

In  particular,  if  x   is  the  turning  point  and  Yy^  ^nd  Y,^  have 
each  been  'connected'  to  the  turning  point,  then  K(e)  is  known 
asjnnptotlcally  with  an  error  term  v/hich  is  of  the  order  of  a 
poijer  of  e.   Consequently,  the  asynptotic  form  of  Y,  has  been 
'continued'  into  the  sector  }Oy,    i.e., 

Y^^Y  in   j/^ 


h 

'VYIC  in   iiJ  ,, 


It  is  illuminating  to  conpare  the  above  discussion  with 
the  follo'-jinc  sketch  of  the  parameter  free  case.  As  x 
approaches  a  non-simple  pole  of  A,  i.e.,  as  x  approaches  an 
lrre;^ular  sinc^ular  point  of  the  equation,  there  exists  a  full 
complement  of  formal  solutions  (Fabry  [1],  Thome  [1], 
Tm'*rittin  [3]).   The  asymptotic  natiore  of  the  formal  solutions 
can  be  shOTn  in  a  rainimvim   sector  (Trjitsinsky  [1])  or  in  a 
maximum  sector  (''asow  [1,  chap.  Ill],  l^iedrichs  [1,  p.  119]) 
or,  for  normal  equations,  in  a  Stokes  sector  (Poincare  [1],  [2], 
3irl-±ioff  [2],   (Of  coTorse,  here  the  Stokes  lines  are  limiting 
directions  rather  than  f inite  lines.)   In  order  to  continue 
as3,Taptotic  forms  into  successive  sectors  it  is  sufficient  to 
'connect'  the  solution  types  of  any  ordinary  point  of  the 
equation  (Turrittin  [1]).   Finally,  one  should  keep  in  mind 
that  as  e  — >  o  one  can  map  a  domain  which  approaches  the 
irregular  singular  point  onto  a  finite  neifjhborhood  of  the 
turning  point.  Thereby  one  solves  a  special  class  of  turning 
point  problems  in  terns  of  a  solvable  irregular  singular  point 
problem.  Actually,  the  irregular  singular  point  is  taken  to 
be  at  X  =  oo  so  that  one  must  stretch  the  variable  of  the 
parameter  problem. 


.r;; :  • 


^'A-r 


III.   The  contour  Intecral  Approach 

A.  The  Triangular  Sj''stem 

If  A  is  a  triangular  r.iatrix,  i.e.,  if  all  the  elements 
of  A  beloi:  the  raain  diagonal  are  zero,  then  the  solution  to 

Y'  =  A(x,e)Y 

can  be  [^iven  explicitly.   More  specifically,  y  x-7ill  be  a 
triangular  matrix  with  exponentials  along  the  main  diagonal. 
e:q30nential  integrals  along  the  first  super  diagonal,  and 
r-tuple  iterated  ex"ponential  integrals  along  the  r+1  super- 
diagonal.   The  exponential  integrals  are  readily  evaluated 
in  a  sector  av;ay  from  the  turning  point  by  what  amounts  to 
the  saddle  point  method.   In  a  similar  manner  one  obtains 
the  connection  to  the  turning  point  and  t'lereby  the  continuation 
into  usccessive  sectors.  Alternately,  the  continuation  problem 
can  be  solved  directly  (Priedrichs  [1,  p.  3-7])« 
Namely,  if  a  path  of  an  exponential  integral  la  from  a^ 

to  X  for  X  in  a  sector    ;<f^n  then  v.hen  ?:  meanders  into  sector 
p^^   one  revrrites  the  path  as  a^  to  ap  and  then  ap  to  x.   The 
integration  from  a^  to  o.^   can  then  be  recognized  as  asymptotic 
to  a  variant  of  the  gamma  fiinction.   The  exponential  integral 
is  thus  continued  asjTnptotically  into  successive  sectors. 
The  singlj"-  iterated  exponential  integral  can  then  be 
asyitiptotically  integrated  with  the  aid  of  the  above  information 
on  the  exponential  integral,  and  so  on. 
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A  point  to  emphasize  is  that  the  continuation  can 
proceed  even  when  a  singularity'-  is  present  at  the  turning 
point.   7or  exaraple. 


has  X  =  0  as  a  turninc  point  but  this  is  an  ordinary  point 
only  if  Y  is  a  non-negative  integer.  Nevertheless, the 
continuation  problem  can  be  solved  in  the  manner  indicated 
above  for  all  y 

An  attractive  feature  of  the  triangular  system  is  the 
ease  with  which  examples  and  counterexamples  can  be  displayed. 
For  example,  the  system 


■yi\'    /I   X 


illustrates  a  case  where  the  canonical  form  of  the  leading 
matrix  alters  abruptly  at  x  =  0  and  yet  a  turning  point  does 
not  occur  at  x  =  0. 

3.   The  Hyperge  one  trie -Type  Eqiiation 

An  equation  vrith  a  parameter  v/hich  is  reduced  to  a 
parameter  free  equation  by  a  simple  stretching  of  the  variable 
must  have  the  follov;ing  form: 


[ 


,   d  \   p   ,   d  \ 


y  =  0 


where  p,  and  Po  are  polynonlals  iritYi   constant  coefficients. 
The  origin  of  the  ::-plane  is  a  tm-'ninc  point  provided  that 

n-  =  decree  of  p,  >  decree  of  ?£  ^  ^2 

(The  characteristic  or  eigenvalue  equation  then  has  n2  roots 
identically  zero  and  n,  -  np  roots  uniformly  distributed  about 
the  turning  point  and  i/hlch  coalesce  to  zero  at_  the  turning 
point. ) 

In  terns  of  the  nev;  (stretched)  variable  t,  where 
et^  =  x^,  the  equation  becomes 

d  X  4.P_  /4.  d 


Pi(t  ^)-^  P2(*  ^)]y  =  0 


Attention  will  now  be  focused  on  solving  the  continuation 
problem  at  the  irregular  singular  point  t  =  co  since  this  is 
sufficient  to  solve  the  continuation  about  the  tuirninc  point 
at  X  =  0.   There  are  noi;  tv;o  decisive  points. 

Firstly'-,  v;hen  this  parameter  free  equations  has  an 
ordinary  point  or  regular  singularity  at  t  =  0,  then  the  non 
poljniomlal  solutions  are  simply  ex;)ressible  in  terms  of 
generalized  hypergeometric  functions.   Secondly,  the  asjTiiptotic 
behavior  of  the  generalized  hypergeometric  functions  is 
sufficiently  v/ell  kno^^^l  to  permit  the  solution  of  the 
continuation  problem  around  t  =  co  .   Further  details  on  the 
second  point  are  in  the  paper  by  Turrittin  [1]  and  only  the 
first  point  will  be  amplified  here. 
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Fuchsian  theory  assures  one  that  t  =  0  is  a  rerular 
singularity  provided 

n,  <  np  +  p    . 

Moreover,  if  none  of  the  roots  of  p,  differ  by  §  multiple  of  p 
then  there  exists  a  full  complement  of  solutions  of  the  form 

7  =  t^I^  a  tP^   .   a^  7^  0 

where  (y  is  a  root  of  the  polj'-nomial  p-,  .  The  coefficients 
satisfy  the  difference  equation 

Pj^  ('^  +  pi+P)a^^^  -  Pj^  (cr+  p/)a^=  0   . 

II oi;  one  need  merely  note  that  this  is  a  two  term  or  binomial 
difference  equation  with  polynomial  coefficients  and  this  in 
turn  can  be  explicitly  solved  in  terms  of  factorials.   Ilore 
precisely,  let  the  polynomials  be  factored  in  the  form 

^1  ^2 

a  .   IT  (  i^  +  p.  )=  a  Y  "TT  (/+  a.  )   . 

Then  ,  „ 

a.  II   II  (k  +  (3  )  =  a  Y^TT   II  (  f  +  a.) 
'^   k=l  i=l     ^     °      k=l  i=l      ^ 

which  can  be  rei-nr-itten  as 

^1  ,^2 

a.  TT(P.  +y.)l/Pil   =ao'^^'n"  (c   +f):/a.:      . 

^i=l    ^  1         °     1=1     ^  ^ 
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Thus  in  the  standard  notation  for  the  generalized  hyper- 
geometric  function 


•^       n^+l  n. 


1,  a^+1;  yt^ 
p,+l 


Finally,  there  is  the  useful  observation  that  the  solutions  of 

[p^  (td/dt)  +  tP  p^  (td/dt)]  y  =  0 

are  contained  among  the  solutions  to 

(td/dt)^{p^  (td/dt)  +  tP  p^  (td/dt) j  y  =  0   . 

Thus  complete  information  with  regard  to  the  former  equation 
yields  partial  information  with  regard  to  the  latter  equation 
and  this  vrith  little  effort. 

C.   Iliscellaneous  Equations 

A  host  of  irregular  singular  point  problems  can  be  and 
have  been  solved  by  the  contour  Integral  approach.   This 
approach  appealed  in  the  preceding  sections  but  uill  novj  be 
depicted  in  a  more  .';eneral  version.  The  approach  requires 
the  solutions  to  have  the  form  of  a  specific  type  of  integral 
(Mellin  integral,  Laplace  integral,  Euler  integral,  Jordan- 
Poclcharamer  integral,  Sorimerfeld  integral,  etc.).   This 
necessitates  firstly  a  restricted  path  of  integration  and 
secondly  an  integrand  which  satisfies  a  transformed  equation 
which  is  either  a  difference  equation  (for  the  Mellin 
integral)  or  a  differential  equation  (for  the 
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others).   Noreover,  a  sufficient  nuriber  of  paths  must  be 
available  so  that  a  complete  set  of  solutions  exist,  each  in 
the  form  of  the  postulated  integral.   Finally,  the  behavior 
of  the  solution  of  the  transformed  equation  must  be  known 
e.'cplicitly  in  the  vicinity  of  each  of  its  singularities. 
Thus  the  customary  transformed  equations  are  the  two- term 
difference  equation  (Scheffe  [1],  Turrittin  [1]),  the  first 
order  differential  equation  (V/asow  [2]),  and  the  differential 
equations  of  the  'classical'  functions  (Linstone  [1], 
'■Jright  [1]). 

!f/hen  the  above  contour  integral  approach  is  not  directly 
applicable  to  an  equation  it  may  become  applicable  after  a  co- 
ordinate transformation.   For  example,  the  Mellin  transform 
of  the  V/eber  equation  ( Ince  [1,  p.  159]) 

y"  +  (n  +  ^  -  ^  x2)y  =  0 

results  in  a  trinomial  difference  equation.  Hov/ever,  the 
change  of  dependent  variables 

V  =  e^   y 
results  in 

v"  -  XV'  +  nv  =  0 

which  is  of  Hypergeometric-type . 


■  %  -X^>'' 


;  ,  f, 
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IV ,   strong  Perturbations 

A.  The  Diagonal  Leading  Matrix 

1.  A  Solution  of  Distinct  Exponential  Type 

This  section  is  concerned  uith  the  vector  equation 

u»  =  (pA  +  B)u 

in  a  finite  neighborhood  of  x  =  0,  79  (x),  vjhere  the  coef- 
ficients are  restricted  by 
Condition  A:   A  is  diagonal 
Condition  B:   B  is  uniforraly  bounded  in  ""^(x)  as  p  — >oo. 

(For  sir.iplicity  B  is  taken  as  independent  of  p  ), 
The  equation  can  be  reviritten  as 

u»  =  (pA  +  D  +  B)u 

vThere 

B  =  D  +  B  . 

The  subsequent  analysis  for  a  particular  solution  requires 
B  to  have  a  particular  term  on  its  riain  diagonal  identically 
zero.   To  simplify  the  discussion  let  D  be  chosen  as  a 
diagonal  raatrix  uith  the  same  diagonal  elements  as  B.   Thus 
all  the  diagonal  elements  of  B  are  novr  zero.  The  solution  of 


Ik 

the  comparison  equation 

G»  =  (pA  +  D)G 

is  simply  the  exponential  of  a  diagonal  matrix.  Moreover 
the  original  equation  is  noivr  equivalent  to  the  integral 
equation 

u(x)  =  g(x)  +  r  G(x)  G'-'-Ca)  B(a)  u(a)  da 

v;here  g  is  a  particular  solution  and  G  is  a  fundamental 
solution  of  the  co;iparison  equation.  For  definiteness  let  g 
be  chosen  as  the  second  column  of  G,  i.e.. 


Si  =  0        i  +  2 


X 

p 
=  exp  '■     ipX^  +   t>22)dt  ,   1  =  2 


where      Xp     is   the   api:>ropriate  eigenvalue  of    A  ,      Let 

X 

v(x)    =  u(x)    exp(-p     I       X2(^)^^) 


ihe  integral  equation  for  v  is 


■  a    -  - 


J-. 
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a 


=  exp  \    bppdt  +  YZ    \    ^i-jV.  da,   1  =  2. 


■    ■'   ' 


or,  in  synbolic  notation 

V  =  CJ  +  Pv 

(Provisionally,  the  presumption  is  that  X.  ^:  Xp  for  i  4  2,  i.e., 
there  is  but  a  single  formal  expansion  of  exponential  type  2) • 
Paths  for  the  integral  equation  are  chosen  in  a  sector  such 
that  the  exponential  terras  in  the  integrand  are  bounded  (by 
unity) ,  as  p  -e»  00  ,  uniformly  for  x  in  the  sector  and  a  along 
the  path.  Thus  in  a  sector  Tjhere  the  asyraptotic  solution  of 
type  2  is  completely  recessive  all  paths  can  be  taken  the  same. 
(For  details  see  Friedrichs  [l,p.l55]).   Of  course,  since  the 
integral  equation  is  of  the  Volterra  type  the  solution  for  all 
finite  p  can  be  exhibited  as  the  infinite  sum  of  the  successive 
iterations  of  the  integral  opero.tor  applied  to  the  inhomo- 
geneous  term,  1.  e.. 


oo 

u  =  2 p'^(  ^  )  • 

r=o 


Nov:  the  key  fact  is  that  the  iterant  solution  of  this  equation 
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forms  an  asyr.iptotlc  series  for  all  x  in  the  sectoi*  including 
X  =  0.   In  the  successive  iterations  it  Is  the  presence  of 
the  exponential  integrals  which  provides  the  rapid  convergence 
of  the  v^,  i  f  2,  according  to  positive  powers  of  e  . 
As  a  typical  exponential  integral  consider 


^o 


exp  p(x^  -  t^)  dt 


in  the  sector  where  the  integrand  is  uniformly  boxmded.  The 
point   |x^I  =  e  ''^   is  called  an  intermediate  point  since 
the  exponential  can  be  regarded  as  highly  dariped  or  very 
oscillatory  only  for   |x|   greater  than   |x. | .  For  this 
portion  of  the  path  one  can  integrate  by  parts  and  obtain 
an  e/  Ixf^^-^  dependence. v/hich  vanishes  at  least  as  rapidly 
as  c  /  -  Along  the  path  from  x  =  0  to  x  =  x.   the  ex- 
ponential is  bounded  while  the  path  length  vanishes  at  least 
as  rapidly  as   e  '   .  Hence 


■ 


exp  p(x^  -  t^-)dt  =  0(e^/^) 


for  all  X  in  the  sector.  Finally^ one  demonstrates  the 
asyraptotic  convergence  of  the  Vp   in  successive  iterant s 
by  modifying  the  norra  of  v.   In  this  modified  norm  C.j  is 
bounded,  the  Vp  converge,  and  the  convergence  of  the 
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remaining  v.   Is  altered  but  not  spoiled.  Roughly,  if 


1/k, 


\  dx  exp  p  I  (X^  -  X2)dt  =  0(e   ^) 
^  'a 


then  one  defines  ^  ~  ^  Min  (l/k.)  and  the  modified  norm 


e 


|v|  =  p"  ri  Iv  I  +  IvJ 

ii2       "^ 


where   jv. (   are  the  maximum  norms  and  concludes  that  in  this 
norm  the  Plcard  iterants  converge  as   e   (Friedrichs  [1,  p. 
19i|-] )  . 

As  a  consequence  of  the  above  key  fact  one  obtains  the 
connection  to  the  turning  point  of  a  particular  solution  by 
successive  Picard  iterants.  (Incidentally,  in  order  that  there 
be  no  exponential  p  dependence  at  x  =  0  it  vjas  necessary 
that  the  eigenvalue  v/as  integrated  from  the  turning;  point 
to  X  ).   One  then  adduces  the  sup  .lemental  fact  that  this 
particular  solution  has,  vjhen  x  is  a  finite  distance  from 
the  turning  point,  precisely  the  desired  asymptotic  form 
proceeding  in  increacing  intef^ral  powers  of   e  =  p"   .   This 
fact  is  clear  from  the  preceding  discussion  of  the  exponential 
integral  where  the  integral  powers  of  e  are  obtained  by 
integration  by  parts.   (Of  course,  it  takes  on  the  order  of 
1/9   iterations  to  obtain  each  successive  term.   Hov;ever, 
the  customary  formal  procedure  already  presents  a  canonical 
form  for  the  asymptotic  expansion) , 


VI' 


C    i'ix 


i      >i 
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2.  A  Solution  of  Non-Distinct  Exponential  Type 

Suppose  one  starts  again  xrith 

u'  =  (pA  +  B)u 

restricted  by  Conditions  A  and  B.   One  seeks  to  connect  to 
the  turning  point  a  solution  of  type  corresponding  to  Ap 
but  novr  ho\Jever  let  X,  =  X^  =  , . .  =  X  i   but  X .  ^  A,   for 

Id  )L  J  '   1 

X  ^  0  and  j  >  ,'/  .   In  order  to  construct  the  appropriate 
comparison  equation  imagine  A  and  B  to  be  partioned  into  four 
parts  by  tvjo  straight  lines,  one  line  just  belovr  rov;  /?.    and 
the  other  line  just  to  the  right  of  column  Jo  ,     B  is  then 
separated  into  the  sum  of  two  parts 


B  =  B-j_  +  B 


A 

where  B  is  identically  zero  in  the  upper  Is  ft  partition 
but  all  other  elements  of  B  are  identically  the  same  as 
those  of  B.  The  conparison  equction  is 


G»  = 


(pA  +  B3_)G 


and  it  is  clear  that  vrith  the  above  partitioning  this  equation 

splits  into  tvro  separate  equations.  Consequently, 

G(x)  G   (a)B  and  B  have  the  same  placement  of  identically 
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zero  elements.  Iloreover,  although  the  upper  left  partition 
of  G(x)   nay  not  be  knovjn  in  closed  form  yet  it  is  vrell 

X 

defined  as  e.xp  p  f  X-,(t)dt  multiplied  by  the  solution  of  a 

o 

parameter  free  linear  ordinary  differential  equation  near  an 

ordinary  point.   Ihe  analysis  proceeds  as  in  the  previous 
section  with  the  modified  norm 


v|  =11  |v.|  +  p®  EZZ  IvJ 

x=l  i=  £+1 
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B.  The  Diagonalizable  Leading  Matrix 

V/hen  the  leading  matrix  is  not  diagonal  it  lo   natural 
to  attempt  to  diagonalize  it  by  a  co-ordinate  transformation, 
If  w  is  the  net-7  dependent  variable  defined  by 

u  =  Pu 

then  the  new  equation  is 


w»  =  (p  p'-'-AP  +  P'-'-BP  -  P'-'-p»)w 


For  the  columns  of  P  one  chooses  independent  eigenvectors 
of  A.   The  eigenvalues  of  A  are  presumed  to  be  completely 
distinct  for  x  4  0  aiid  partially  coalescent  at  x  =  0. 
The  matrix  p"  AP   has  then  the  eigenvalues  of  A  along  the 
main  diagonal  and  zeros  elseuhere.   These  eigenvp.lues  are 
certainly  continuous  at  x  =  0.  Moreover^  if  cj)  is  a 
coliimn  of  P  corresponding  to  the  eigenvalue   X  then  the 
components  of  (j)  are  proportional  to  the  cofactors  of  some 
row  in  A-XI.  Hence  the  elements  of  P  can  certainly  be 
chosen  to  be  continuous  at  x  =  0.   For  the  moment  let  us 
presume  that  P"   is  continuous  at  x  =  0.   Then  tlie  new 
equation  can  be  written  as 

w»  =  (pD  +  B)w 
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v;here  D  is  diagonal  vjith  Lipschitz  continuous  elements  and 
XxThere  the  elements  of  B  are  integrable ,   The  argument  now 
proceeds  as  given  in  the  preceding  section  nith  obvious 
alterations.  For  example,  the  asymptotic  estimate  of  each 
exponential  integral  term  novi   depends  upon  the  degree  of 
integrability  of  the  integrand,  i.e.,  for  q  >  -1 

^o  1+q 


\      t^   e::p  p(x^^  -  t^)dt  =  0  (e""^) 


for  all  X,  including  x  =  0,     in  an  appropriate  sector.  Also 
the   9  to  be  employed  for  the  modified  norm  is 


e  =  i  Min(  ^   )  . 


Let  us  novT  briefly  return  to  the  restrictions  to  be 
imposed  upon  A(x)  so  that  both  P(x)   and  P   (x)   are 
continuous  at  x  =  0  uhere   P"-^AP  is  diagonal.   Obviously, 
it  is  necessary  that  A(0)   have  simple  elementary  divisors. 
More  precisely,  it  is  necessary  that  the  derivative  of  lovrest 
order  at  x  =  0,  say  A^^'(O),  nhich  is  not  a  scalar  matrix 
have  simple  elementary  divisors.   (This  form  of  the  necessary 
condition  is  apparent  from  the  trivial  remark  that  P  will 
diagonalize  any  scalar  matrix) .   On  the  other  hand,  it  is 
sufficient  that  the  above  defined  A^'^'(O)   have  distinct 
eigenvalues  (since  the  columns  of  P,  i.e.,  the  eigenvalues 
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of  A,  will  then  certainly  remain  independent).  Alter- 
nately, it  is  sufficient  that  A(x)   be  normal  (for  then 
the  normalized  eigenvectors  are  orthogonal  rjid  hence  cannot 
coalesce) . 
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V .  VJeak  Perttrpbati ons 

The  vreak  perturbation  dispenses  vjith  the  estimates 
for  the  exponential  integrals,   it  is  enough  that  a  fundamental 
solution  of  the  comparison  equation  be  bounded  as  p  — >  oo 
uniformly  for  x  on  a  closed  line  from  the  turning  point. 
The  rapid  convergence  of  the  iterants  of  the  integral  equa- 
tions is  then  provided  by  the  portion  of  the  integral  kernel 
which  corresponds  to  the  vreak  perturbation.  The  t;eaker  the 
perturbation,  the  more  rapid  the  convergence  for  large  p  • 
To  complete  the  picture  one  should  delimit  the  class  of 
comparison  equations  vrhich  satisfy  the  boundedness  condition. 
A  possible  approach  to  the  problem  is  the  following.   The 
classical  demonstration  of  the  asyraptotic  nature  of  the  formal 
solutions  first  sets  up  a  Volterra  integral  equation  whose 
inhomogeneous  terra  is  a  truncated  formal  solution.  This 
integral  equation  is  then  solved  by  iteration  for  the  variable 
limit  of  the  integral  within  a  sector  but  alvrays  s  finite 
distance  from  the  turning  point,   Mowever,  one  can  recognize 
the  successive  iterants  to  be  the  ubiquitous  iterated  ex- 
ponential integrals.  Furthermore,  the  inhomogeneous  term 
admits  a  modified  norm  and  finally  the  appropriate  element 
of  the  kernel  is  asj'nptotically  small.   Hence  by  an  application 
of  the  strong  perturbation  scheme  one  extends  the  far  zone 
(vjhere  the  formal  expansions  are  asyi^iptotic)  doi.Ti  to  an 


Zk 


Inflnltesmal  distance  (n.s   p  — ^  oo)  from  the  turning  point. 
But  in  the  near  zone  a  sinple  stretching  of  the  variable 
results  in  an  equation  vjith  coefficients  xonifonnly  bounded 
as  p  — >  00  •  Hence  in  the  near  zone  the  solution  is  a  bounded 
function  of  the  initial  values,  ''hen  the  near  and  far  zones 
raeet  at  an  intermediate  point  then  a  bounded  comparison 
solution  can  be  constructed  . 
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VI,   The  Decoupling  or  Reduc tlon  A pproa c h 

The  decoupling  approach  atteinpts  to  regard  a  system 
as  a  relatively  v;eak  coupling  of  coriiponents  where  the  behavior 
of  each  component  is  either  Icnovm  or  readily  obtainable.  From 
this  one  can  usually  deduce  the  behavior  of  the  coupled 
system.  For  example,  suppose 


G 


11  =  P^ll^il 


and 


'22 


pA^^G 


22"  22 


where  A,,  ,  App   are  square  matrices  and  G-, -,  ,  Gpp  are 
fundamental  matrix  solutions.   Then  a  solution  to 


*n 


=  <p 


0 

+ 

22 

B 

B 


12 


21   ^ 


Is 


^1  =  Si  +  ,3  ^2.1   ^12  ^ 


U2  =  g2  +  j  G22  B21  U2 
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If  now  g  =  g^'Tjg^   is  a  colui-in  of  the  coi'iposite 

G  —  G-j_2  (t) ^22  ^^"^  ^^   ^^®  appropriate  as^i-.iptotic  behavior 
of  G  is  knoxm  then  the  problem  is  ideally  set  up  for  the 
strong  perturbation  formalism. 

As  a  second  and  final  exanple,  let  us  consider  the 
decoupling  due  to  non-coalescing  eigenvalues  (Langer  [2]). 
The  equation  is  as  before 

u»  =  (pA  +  B)u 

where  A (x)   has  distinct  eigenvalues  for  x  4  0  but  now 
one  eigenvalue,  say  X^ ,  remains  distinct  at  the  turning 
point  X  =  0.  For  siraplicity  X^  can  be  taken  as  identically 
zero.   Let   [^(i\]   ^le  a  truncated  formal  solution  corresponding 
to  A^  .  The  key  fact  is  that   t^n )  ^   ^^  analytic  in  a  full 
neighborhood  of  the  turning  point.   This  fact  can  be  seen  as 
follov;s.  Since  the  eigenvalue  remains  distinct  it  certainly 
is  analytic  and  consequently  the  eigenvector  can  be  chosen 
to  have  analytic  components.   One  augments  this  vector  by 
additional  analytic  vectors  to  form  a  complete  unitary  set. 
Then  one  transforms   u  by  this  vector  set  regarded  as  a 
matrix.   The  leading  matrix  of  the  transfon-ied  equation  vrill 
have  only  zeros  in  the  first  coluron.  "hen  a  formal  solution 
of  type  1  is  substituted  into  this  transformed  equation  the 
equations  which  govern  the  successi/e  tem-.is  reveal  these 
successive  terms  to  be  analytic. 


K-  "-'i'^^is'^  V. 
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since   [u. ^.]   is  an  approximate  solution  even  at  the 
turning  point  hence  it  is  not  surprisin.^  that  one  can  reduce 
the  order  of  the  original  equation  in  the  standard  manner 
(HasoTv  [1,  p.  28]).   Namely,  define  P  from  the  equation 
P[u,,,]  =  (o-.)   and  then  chanf;e  the  dependent  variable  to  w 
by  the  transf ornation 

Pu  =  w 

The  equation  for  w  is  then 

,-1  .  ^.-,-1" 


w 


'  =  {_P(pA  + 3)P"  +P'P"  }  w  =:  nr   . 


But  from  the  definition  of  P,  vj  has  a  solution  asymptotic 
to  (6^.)«   Consequently,  the  first  column  of  P  is  asymptotic 
to  zero  so  that  the  last   n-1  components  of  w  decouple 
asymptotically.   Of  co\irse,as  compared  nith  the  full  system, 
the  decoupled  system  has  less  transition  lines  and  hence 
larger  maximum  sectors.  Suppose  then  that  a  cor.iplete  set  of 
vector  solutions  of  the  decoupled  equation  is  connected  to 
the  turning-  point  in  this  larger  sector.   These  vectors 
become  solutions  of  the  full  equation  v;hen  the  appropriate 
first  component  of  i:     is  appended.  These  first  components 
are  readily  determined.  However,  they  are  then  connected  to 
the  turning  point  in  a  maximum  sector  determined  by  the  full 
equation. 
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